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Abstract 

We consider the Ricatti equation in the context of population dy- 
namics, quantum scattering and a more general context. We examine 
some exactly solvable cases of real life interest. 

1 Introduction 

The Riccati Equation [TJ |3 El 

y' = p(t)+q(t)y + r(t)y 2 (1) 

appears in several branches of applicable mathematics, for example popu- 
lation dynamics and mathematical physics, as in quantum scattering.lt is 
known that it exhibits chaotic behavior jl].It is also well known that, it 
reduces to a second order linear differential equation by the substitution, 

y = -, 

when we get, 

u" + g{t)u' + f{t)u = (2) 

(1) and (2) have been studied in detail for a long time [SI El E|. 
It may be mentioned that if two solutions of the Riccati equation are equal at 
a point, then they coincide. The reason lies in the continuity of the logarithmic 
derivative — given in the substitution that lead from (1) to (2). 
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We will now consider exact solutions of the Riccati equation in two different 
contexts, from the field of population dynamics and the field of quantum 
scattering. 



2 Some exact solutions 

2.1 A Problem In Population Dynamics 

We now examine coupled Riccati type equations of relevance to a population 
model. 



Gause's Model 

Consider two-species populations occurring together, and assume that the 
growth of each is inhibited by members ,both of its own and of the other 
species. Denoting the number of individuals in species 1 as Ni and species 2 
as N 2 ,we have the Gause's competition equations: |8] 

1 dN x 

n - aiJVx - a 2 N 2 (3) 



Ni dt 

1 dN 2 
W 2 ~dT 



r 2 - a 3 Nt - a±N 2 (4) 



where ri,r 2 ,ai,a 2 are defined below. Thus we are assuming that the per capita 
growth of each population at an instant is a linear function of the sizes of 
the two competing populations at that instant. Each population would grow 
logistically if it were alone with logistic parameters r\ and a± for species 1 
and r 2 and a 2 for species 2 . 

In general the simultaneous differential equations cannot be explicitly solved. We 
now consider a particular set of circumstances in which they can be solved. Thus 
for (3) and (4) we specialize to, with a more convenient notation , 



dy 
dt 

dx 
~dt 



xy 2 (5) 
yx 2 (6) 
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From (5) and (6) we get, 
dy 1 — xy 2 



(7) 



dx 1 + yx 2 
Integrating (7) we have 

y-x = -x 2 y 2 + C (8) 

From (8) when y = 0,C = — Xo so that y < x 
Reverting back to N\ and N2, this is 

N x - N 2 = -N 2 2 N 1 2 + C (9) 

From (9) when N x = 0,C = -N so that iVi < N 2 . 

The Figure below illustrates a particular case of the above solution with the 
two populations N% and N 2 along Y and X axes. 



2.2 Problems from Scattering Theory 

By reversing the method used above ,the quantum mechanical radial Schrodinger 
equation 

u" - H(r)u = (10) 

can be reduced to a Riccati equation 

v' + v 2 = H{r) (11) 

by the substitution v = — . The form (10)or (11) is used in phase shift analy- 
sis for example in Calogero's variable phase approach j2j.lt is also possible to 
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use the form (11) for building up an iterative procedure [TU] .We would now- 
like to point out that (11) can be used in a different context for providing 
exact solutions for specified classes of the potential function f(r). Let us write 
in (11) 

1 

r = i - a/ 

Whence we get 

n npJ 9dr 

H= l = _J& (12) 

/ J eJ 9dr + C 

l.We put g=a in (12) 

h « 2 



1 + cae~ ar 
The potential H is now of the form 

H ° 



A + Be ar 

This is the well-known Wood-Saxon potential. The solution u is given by 

u = C 1 (e Qr + C 2 a) 
2. Putting g=r in (12), the potential H is given by 



H 



re r 2 /2 



/ e r2 / 2 dr + c\ 

This is a modified Gaussian potential and the solution u is given by 

u = K J e~ r2 dr + C 3 

The Graph is shown in the figure below and is seen to fall very steeply 
indicative of a confined state or particle. 



4 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



-20 -10 ^ TO 20~ 

3. Putting g=l/r in (11), the potential H is given by 

H 2 



r z + c 4 

This is a shifted inverse square potential and the solution u is given by 
u = K'{r 2 + C 4 ) 

3 Other Exact and Asymptotic Solutions 

We now consider some asymptotic solutions of (11), which we write as, 

u' = u 2 + f (13) 

Let 

/ fdr = 9, 

f and g being bounded functions. So (13) can be written as , 

z' = (z + gf > 0, z = u - g, (14) 

which shows that z is an increasing function of r. Suppose z is unbounded. So 
for large r,we should have from (14) 

2 

z m z 
whence 



which — > as r — > oo. 

This is a contradiction. Therefore z is bounded and so also u ,that is u — > M 
as r — ► oo 

Therefore u' ~ (asymptotically) 
Therefore for large r,(13) becomes 

u 2 = —f, whence 

u = ±{-ff/ 2 (15) 
By way of verification of (15), let us consider (13) with, 

So we expect that for large r, u ~ ± - by (15). Let us put 

u = ^ (16) 
r 

in (13). So 

/? f3 2 a2 



fy*^i fy*^i /y*^! 



or 



a -1± Vl + 4a 2 

/? = w ±a if a > 1 

A, 

That is 

a 

K = ± — 

r 

(everywhere, and so also for large r). Using this example,with transfor- 
mations of the independent variable,we can generate similar solutions. For 
example, if we substitute for r, t — t(r), (13) becomes 

t'u = u 2 + f, 

where the dot denotes the derivative with respect to t.The choice t' = f(r), 
leads to a similar equation, and one can verify that for 

f( r ) = %- 2 / 3 - a 2 r 2 '\ u = ar 1 ^ 

is a solution. More generally as can be easily verified ar n is a solution for 
f(r) = nar n ~ 1 — a 2 r 2n and so on. However in these examples, neither u nor 
u' are asymptotically bounded. 

Finally it maybe observed that if one solution of the Riccati equation (13) is 
known, then others could be derived therefrom[3 j . 
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